Metastability and Coherence of Repulsive Polarons in a Strongly Interacting Fermi Mixture 
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Ultracold Fermi gases with tuneable interactions repre- 
sent a unique test bed to explore the many-body physics 
of strongly interacting quantum systems [1-4]. In the past 
decade, experiments have investigated a wealth of intrigu- 
ing phenomena, and precise measurements of ground- 
state properties have provided exquisite benchmarks for 
the development of elaborate theoretical descriptions. 
Metastable states in Fermi gases with strong repulsive in- 
teractions [5-11] represent an exciting new frontier in the 
field. The realization of such systems constitutes a major 
challenge since a strong repulsive interaction in an atomic 
quantum gas implies the existence of a weakly bound 
molecular state, which makes the system intrinsically un- 
stable against decay. Here, we exploit radio-frequency 
spectroscopy to measure the complete excitation spectrum 
of fermionic 40 K impurities resonantly interacting with a 
Fermi sea of 6 Li atoms. In particular, we show that a well- 
defined quasiparticle exists for strongly repulsive interac- 
tions. For this "repulsive polaron" [9, 12, 13] we measure 
its energy and its lifetime against decay. We also probe 
its coherence properties by measuring the quasiparticle 
residue. The results are well described by a theoretical 
approach that takes into account the finite effective range 
of the interaction in our system. We find that a non-zero 
range of the order of the interparticle spacing results in 
a substantial lifetime increase. This major benefit for the 
stability of the repulsive branch opens up new perspec- 
tives for investigating novel phenomena in metastable, re- 
pulsively interacting fermion systems. 

Landau's theory of a Fermi liquid [14] and the underlying 
concept of quasiparticles lay at the heart of our understand- 
ing of interacting Fermi systems over a wide range of energy 
scales, including liquid 3 He, electrons in metals, atomic nu- 
clei, and the quark-gluon plasma. In the field of ultracold 
Fermi gases, the normal (non-superfluid) phase of a strongly 
interacting system can be interpreted in terms of a Fermi liq- 
uid [15-18]. In the population-imbalanced case, quasiparti- 
cles coined Fermi polarons are the essential building blocks 
and have been studied in detail experimentally [16] for attrac- 
tive interactions. Recent theoretical work [9, 12, 13] has sug- 
gested a novel quasiparticle associated with repulsive interac- 
tions. The properties of this repulsive polaron are of funda- 
mental importance for the prospects of repulsive many-body 
states. A crucial question for the feasibility of future experi- 
ments is the stability against decay into molecular excitations 
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FIG. 1 : Energy spectrum of the impurity in the Fermi sea. For 

the two polaronic branches, the energies E+ (red line) and E- (green 
line) are plotted as a function of the interaction parameter — \/{Kpd) . 
The shaded area between the dashed lines representing E m and 
E m — £f ( see t ex shows the continuum of molecular excitations. 
The vertical lines at \/{Kpd) = ±1 indicate the width of the strongly 
interacting regime. The inset illustrates our rf spectroscopic scheme 
where the impurity is transferred from a noninteracting spin state |0) 
to the interacting state 1 1). 

[11, 12, 19]. Indeed, whenever a strongly repulsive interaction 
is realized by means of a Feshbach resonance [20], a weakly 
bound molecular state is present into which the system may 
rapidly decay. 

Our system consists of impurities of fermionic 40 K atoms 
immersed in a large Fermi sea of 6 Li atoms, which is char- 
acterised by a Fermi energy Ep = h X 37 kHz and a temper- 
ature T = 0.16£f/A:£ (see Methods), with h and kg denoting 
Planck's and Boltzmann's constants. In a particular combina- 
tion of spin states [21], the 6 Li - 40 K mixture features a Fesh- 
bach resonance centered atfio = 154.719(2)G. The resonance 
allows to widely tune the s-wave interaction, parametrised 
by the scattering length a, via a magnetic field B. The in- 
teraction strength is described by the dimensionless parame- 
ter — l/(Kpa), where K F = h^ l y / 2mu£F = 1 /(2850ao) is the 
Fermi wave number; here h = h/2n, «o is Bohr's radius, and 
mu is the mass of a 6 Li atom. Near the resonance center, 
the linear approximation — l/(jC/?a) = (B — Bo) /20mG holds. 
The momentum dependence of the interaction is characterised 
by the effective range, which we express in terms of the pa- 
rameter [22] R* = 2700ao (see Supplementary Information). 

Figure 1 illustrates the basic physics of our impurity prob- 
lem in the T = limit, showing the energies of different states 
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as a function of the interaction parameter. The situation is 
generic for any impurity in a Fermi sea, but quantitative de- 
tails depend on both the mass ratio and the particular width 
of the Feshbach resonance. The theoretical curves are based 
on an extension of an approach presented in Refs. [12, 23] 
to our case of a relatively narrow Feshbach resonance with 
KfR* = 0.95 and thus a considerable effective range of the 
interaction (see Supplementary Information). 

The spectrum exhibits two quasiparticle branches, which do 
not adiabatically connect when the resonance is crossed, and 
a molecule-hole continuum (MHC). The interaction-induced 
energy shifts of the two branches (E + > and £_ < 0) are 
generally described in a many-body picture by dressing the 
impurities with particle-hole excitations. Far away from res- 
onance this simplifies to a mean-field shift proportional to a. 
The lower branch £_ of the system (green line) corresponds 
to the attractive polaron, which has recently received a great 
deal of attention theoretically [4, 15, 23-25] as well as ex- 
perimentally [16, 17, 26]. This polaronic branch remains the 
ground state of the system until a critical interaction strength 
is reached, where the system energetically prefers to form a 
bosonic molecule by binding the impurity to an atom taken 
from the Fermi sea [23, 25, 27]. The continuum arises from 
the fact that a majority atom with an energy between and £f 
can be removed from the Fermi sea to form the molecule. The 
MHC thus exists in an energy range between E m and E m — 
(dashed lines in Fig. 1), where E m represents the energy of 
a dressed molecule including the binding energy of a bare 
molecule in vacuum and a positive interaction shift. The at- 
tractive polaron can decay into a molecular excitation if this 
channel opens up energetically (£_ > E m — £f). 

The upper branch (red line in Fig. 1) corresponds to the re- 
pulsive polaron [9, 12, 13] with an energy E + > 0. Approach- 
ing the resonance from the a > side, E + gradually increases 
and reaches a sizeable fraction of £f. However, the polaronic 
state becomes increasingly unstable as it decays to the lower 
lying states (attractive polaron and MHC). Close to the res- 
onance center, the polaronic state becomes ill-defined as the 
decay rate approaches E + /%. 

To investigate the excitation spectrum of the impurities, we 
employ radio-frequency (rf) spectroscopy [28-30]. We ini- 
tially prepare the 40 K atoms in a non-interacting spin state |0) 
and then, with a variable frequency Vrf , drive rf transitions into 
the resonantly interacting state |1). Our signal is the fraction 
of atoms transferred, measured as a function of the rf detuning 
Vrf — Vo with respect to the unperturbed transition frequency 
Vo between the two spin states. This excitation scheme pro- 
vides access to the full energy spectrum of the system. In par- 
ticular, it allows us to probe the metastable repulsive polaron 
as well as all states in the MHC. We furthermore take advan- 
tage of the coherence of the excitation process by driving Rabi 
oscillations. As an important practical advantage, this enables 
very fast and efficient transfer of population into a short-lived 
quasiparticle state by application of 7T -pulses. Moreover, we 
will show that measurements of the Rabi frequency directly 
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FIG. 2: Spectral response of 40 K impurities in a 6 Li Fermi sea. 

The false-colour plots show the fraction of 40 K atoms transferred 
from the non-interacting spin state |0) into the interacting state |1) 
for different values of the rf detuning parameter A = ^(Vrf — Vo) and 
for variable interaction strength —l/(Kpd). The panels a and b refer 
to low and high rf power. For comparison, the lines correspond to 
the theoretical predictions for E + , E,„, and E,„ — £p as shown in 
Fig. 1. 



reveal quasiparticle properties. 

In Fig. 2 we show false-colour plots of our signal, detected 
for different values of the detuning parameter A = /j(Vrf — Vo) 
and for variable interaction strength — 1 / ( Kf a) . Figure 2a dis- 
plays a set of measurements that was optimised for signal and 
spectral resolution of the polaronic excitations by using mod- 
erate rf power (see Methods). The spectrum in Fig. 2b was 
optimised for detection of the molecular excitations. Here 
a much higher rf power had to be employed because of the 
reduced Franck-Condon wavefunction overlap. For the pola- 
ronic branches the high rf power leads to a highly nonlinear 
saturation behaviour. 

Our data clearly show both polaronic branches, with their 
measured energies being in excellent agreement with theory. 
The attractive polaron is found to disappear in the strongly in- 
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FIG. 3: Decay rate of the repulsive polaron. The data points dis- 
play the measured decay rates T as extracted by exponential fits from 
decay curves; the error bars indicate the fit uncertainties. Sample de- 
cay curves are shown in the inset. The solid lines represent theoreti- 
cal calculations of the two-body decay (blue line) and the three-body 
decay (red line) into the attractive polaron or the MHC, respectively. 



teracting regime. This behaviour, which is different from the 
one observed for 6 Li spin mixtures [16], is consistent with the 
crossing of E- and E m — Ef at —1/ (jQ?a) « +0.6 predicted 
for our system. In contrast, the repulsive polaron extends far 
into the strongly interacting regime. A sharp peak is observed 
in the spectrum with decreasing signal strength until it finally 
fades out very close to resonance at — l/(K)?a) ~ —0.3 (see 
Supplementary Information). The low rf power only produces 
a weak signal for the MHC, whereas the high rf power clearly 
unveils the presence of the MHC. Further on the a > side 
of the resonance, the molecular signal gets weaker because 
of the decreasing Franck-Condon overlap, and outside of the 
strongly interacting regime the situation corresponds to the rf 
association of bare molecules (see Supplementary Informa- 
tion). 

To investigate the decay of the repulsive branch, we apply 
an rf pulse sequence to selectively convert repulsive polarons 
back into non-interacting impurities after a variable hold time 
(see Methods). The back conversion sensitively depends on 
the rf resonance condition and thus allows us to discriminate 
40 K atoms in the polaronic state against those ones forming 
molecules. Figure 3 presents the experimental results. The 
inset shows three sample curves, taken for different values 
of the interaction parameter. The main panel displays the 
values extracted for the decay rate F from the decay curves 
by simple exponential fits. The data reveal a pronounced 
increase of decay as the resonance is approached, which is 
in good agreement with theoretical model calculations [12] 
(see Supplementary Information). The decay populates the 
MHC and may happen in a two-step process where the re- 
pulsive polaron first decays via a two-body process into an 
attractive polaron (blue line) and then decays into a molecular 
excitation. Alternatively, a three-body process may directly 
lead into the MHC (red line). Very close to the resonance, 
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FIG. 4: Rabi oscillations and the quasiparticle residue, a, Sample 
Rabi oscillations (magenta and green data points for — 1/(kj?o) = 
— 1.25 and —0.5, respectively) with harmonic oscillator fits (solid 
lines) demonstrate the two effects of the interaction with the Fermi 
sea: damping and a reduction of the Rabi frequency. The black curve 
is a reference curve recorded without 6 Li . In b and c, the data points 
show the damping rates y and the normalized Rabi frequencies Q,/Qq 
as measured for two different values of the rf power; the blue squares 
and red dots refer to Q.q = 2tz x 6.5 kHz and 12.6 kHz, respectively. 
The error bars indicate the fit uncertainties. The solid lines represent 
the theoretical behavior of yZ for the repulsive and the attractive 
polaron. 



for —1/(Kfo) = -0.25, we find TiT/ef ~ 0.01, which corre- 
sponds to a l/e lifetime of about 400 jJ.s. Relating this decay 
rate to the corresponding energy shift E + = Q.3QEf, we obtain 
hT/E + w 0.03 <C 1, which demonstrates that the repulsive po- 
laron exists as a well resolved, metastable quasiparticle even 
deep in the strongly interacting regime. 

The lifetime observed for the repulsive branch appears to 
be remarkably long, when compared to recent experiments on 
6 Li spin mixtures [11]. The latter system is a mass-balanced 
one and it features a broad Feshbach resonance with a negli- 
gible effective range. Our theoretical approach allows to an- 
swer in a general way the question how mass imbalance and 
the width of the resonance influence the lifetime. We find that, 
while the mass imbalance does only play a minor role [12], the 
dominant effect results from the finite effective range, which 
is associated with the narrow character of the Feshbach res- 
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onance that we exploit. Comparing our situation with a hy- 
pothetical system with a broad Feshbach resonance, and thus 
with a zero-range interaction, we find that in the strongly in- 
teracting regime the same amount of energy can be obtained 
with an almost ten times increased lifetime (see Supplemen- 
tary Information). 

Besides energy and lifetime, the polaron is characterized 
by its effective mass m* and its quasiparticle residue Z. The 
difference between effective and bare mass [12] does not pro- 
duce any significant features in our rf spectra. The residue Z 
(0 < Z < 1) quantifies how much of the non-interacting par- 
ticle is contained in the polaron's wavefunction, which can 
be written as \fX | 1 ) plus terms describing excitations in the 
Fermi sea. The pre-factor \[Z directly manifests itself in the 
Rabi frequency £2 that describes the coherent rf coupling be- 
tween the noninteracting and the polaronic state (see Supple- 
mentary Information). 

Figure 4 presents the experimental data on Rabi oscillations 
for variable interaction strength. The sample curves in Fig. 4a 
demonstrate both the interaction-induced change in the fre- 
quency and a damping effect. We apply a simple harmonic 
oscillator model (including a small increasing background) to 
analyse the curves, which yields the damping rate y and the 
frequency £2. The damping strongly increases close to the res- 
onance center, but does not show any significant dependence 
on £2o, see Fig. 4b. It is interesting to note that the population 
decay rates T measured for the repulsive branch (Fig. 3) stay 
well below the values of y, which points to collision-induced 
decoherence as the main damping mechanism. 

Figure 4c displays the measured values for the Rabi fre- 
quency £2, normalized to the unperturbed value £2o. The 
interaction-induced reduction of £2/£2o is found to be indepen- 
dent of the particular value of £2q (comparison of blue squares 
and red dots; see also Supplementary Information). The solid 
lines show \[X as calculated within our theoretical approach 
for both the repulsive and the attractive polaron. The com- 
parison with the experimental data demonstrates a remarkable 
agreement with the relation \fZ = D./D.Q. Our results there- 
fore suggest measurements of the Rabi frequency as a precise 
and robust method to determine the quasiparticle residue Z, 
and thus provides a powerful alternative to methods based on 
the detection of the narrow quasiparticle peak in the spectral 
response [16, 23]. 

In conclusion, we have realized an ultracold model system 
of 40 K and 6 Li atoms to investigate the quasiparticle behavior 
of heavy impurities resonantly interacting with a Fermi sea 
of light particles. Our spectroscopic approach has confirmed 
the existence of the predicted repulsive branch [9, 12, 13] 
and has demonstrated that the repulsive polaron can exist as 
a well-defined quasiparticle even deep in the strongly inter- 
acting regime. The long lifetime of the repulsive polaron in 
our system, which we ascribe to the finite effective range of 
the interparticle interaction, may be a key factor to overcome 
the problem of decay into molecular excitations [11, 19] in 
the experimental investigation of metastable many -body states 
that rely on repulsive interactions. In particular, the creation 



of states with two fermionic components phase-separating on 
microscopic or macroscopic scales [5-8, 10, 1 1] appears to be 
an intriguing near-future prospect. 
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METHODS 

Experimental conditions. Our system consists of 2 x 10 4 
40 K atoms and 3.5 x 10 5 6 Li atoms confined in an optical 
dipole trap. The trap is realized with two crossed beams de- 
rived from a 1064 nm single-mode laser source. The measured 
trap frequencies for Li (K) are v, = 690 Hz (425 Hz) radially 
and V- = 86 Hz (52 Hz) axially; this corresponds to a cigar- 
shaped sample with an aspect ratio of about 8. The preparation 
procedure is described in detail in Ref. [31]. The Fermi ener- 
gies, according to the common definition Ef = h %/6Nv}v z 
for harmonic traps, are Ep 1 = h X 44kHz = kg x 2.1 jj,K and 
Ef = h X 10.4kHz = k B x 500nK. At a temperature T ps 
290 nK the 6 Li component forms a deeply degenerate Fermi 
sea (ksT /Ep 1 pa 0. 14) while the 40 K component is moderately 
degenerate (k B T/Ef ps 0.6). 

Effective Fermi energy. The 40 K atoms experience a 
nearly homogeneous 6 Li environment. This is because the 
optical trapping potential for 40 K is about 2 times deeper than 
for 6 Li and the 40 K cloud is confined in the center of the much 
larger 6 Li Fermi sea [32]. This allows us to describe the sys- 
tem in terms of the effective Fermi energy Ef , defined as the 
mean Fermi energy experienced by the 40 K atoms. We find 
Ef = hx 37 kHz, with two effects contributing to the fact that 
this value is about 15% below Ep, The finite temperature re- 
duces the Li density in the trap center, leading to a peak local 
Fermi energy of h x 40kHz. Moreover, the 40 K atoms sample 
a small region around the trap center, where the density and 
local Fermi energy are somewhat lower. The distribution of 
Fermi energies experienced by the 40 K cloud, i.e. the residual 
inhomogeneity of our system, can be quantified by a standard 
deviation of h x 1 .9kHz. 

Concentration. The mean impurity concentration (mean 
density ratio «K/ n Li) is about 0.4, if one considers the pop- 
ulation of K atoms in both spin states. This may be a priori 
too large to justify the interpretation of our data in terms of 
the low-concentration limit. We find that this interpretation 
is nevertheless valid, as we take advantage of several facts. 
Under strongly interacting conditions only a fraction of the 
K atoms is transferred into spin state |1) (see Fig. 2), which 
reduces the concentration of interacting impurities. A recent 
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quantum Monte Carlo calculation of the equation of state of 
a zero-temperature 6 Li - 40 K Fermi-Fermi mixture [33] fur- 
ther supports our interpretation in the low-concentration limit: 
The strongest interaction in the mass-imbalanced system is 
expected when one has about 4 times more 40 K atoms than 
6 Li atoms, and for concentrations up to a value of 1 the inter- 
action energy per 40 K atom is expected to remain essentially 
constant. To support our basic assumption with experimen- 
tal data, we also took rf spectra for variable numbers of 40 K 
atoms, confirming that in the relevant parameter range finite 
concentration effects remained negligibly small. 

Interaction control via Feshbach resonance. The Fesh- 
bach resonance used for interaction tuning is discussed in de- 
tail in Refs. [21, 32]. It is present for 6 Li in the lowest spin 
state and 40 K in the third-to-lowest spin state. The latter rep- 
resents our interacting state |1); the corresponding quantum 
numbers are F = 9/2 for the hyperfine and mp = — 5/2 for the 
magnetic sub-state. The neighboring state with mp = —7/2 
serves as state |0); here the interspecies scattering length 
(about +65 ao with «o being Bohr's radius) is so small that it 
can be neglected to a good approximation. The tunable scat- 
tering length for state 1 1) in the Fermi sea is well described 
by the standard formula a = «b g ( 1 — AB /(B — Bq)) with «bg = 
63.0a , AB = 880mG, and B = 154.719(2) G. Note that the 
value given for Bq refers to the particular optical trap used in 
the experiments, as it includes a small shift induced by the 
trapping light. The value therefore somewhat deviates from 
the one given in Refs. [21, 32]. In free space, without the 
light shift, the resonance center is located at 154.698(5) G. 
The character of the resonance is closed-channel dominated 
[20]. Following the definition [22] of a range parameter 
R* = /i 2 /(2m,-flbg 8[l AB), with m r = mumx/ '(mu + m^) be- 
ing the reduced mass and 8fi the differential magnetic mo- 
ment, the resonance is characterized by R* = 2700 ao- This 
value accidentally lies very close to l/Kp = 2850 ao, which 
also means that the strongly interacting regime roughly cor- 
responds to the universal range of the resonance. Our sys- 
tem therefore represents an intermediate case (KpR* = 0.95), 
where the behavior is near universal, but with significant ef- 
fects arising from closed-channel contributions. 

Details on Rf pulses. For taking the data of Fig. 2 we 
used Blackman pulses [34] to avoid side lobes in the spec- 
trum. For the upper panel, the pulses were 1 ms long (spec- 
tral width 0.7 kHz ~ 0.02 Sp /h) and the rf power was adjusted 
such that K -pulses would be realized in the absence of interac- 
tions with the Fermi sea. For the data in the lower panel, the 
rf power was increased by a factor of 100 and the pulse du- 
ration was set to 0.5 ms. This resulted in pulses with an area 
of 5n without the Fermi sea. For the lifetime measurements 
in Fig. 3, we used a sequence of 3 Blackman pulses. The first 
pulse (duration between 150/is and 500 /is) was set to drive 
the non-interacting impurity from spin state |0) (mp = —7/2) 
into state |1) (mp = —5/2); here the frequency was carefully 
set to resonantly create repulsive polarons and the pulse area 
was set to fulfill the n -pulse condition. The second pulse 
was a short (60 /is) cleaning pulse, which removed the pop- 



ulation remaining in |0) by transfer to another, empty spin 
state (mp = —9/2). The third pulse had the same parameters 
as the first one and resonantly back-transferred the population 
from the polaronic state in |1) to the non-interacting state |0), 
where it was finally measured by spin-state selective absorp- 
tion imaging. The measurements of Rabi oscillations in Fig. 4 
were performed with simple square pulses. 



SUPPLEMENTARY INFORMATION 
1. Theoretical framework 

The theoretical results presented in the main text and in this 
Supplementary Information are obtained from a model that 
describes the behaviour of a single impurity embedded in a 
Fermi sea with tuneable s-wave interaction near a Feshbach 
resonance with arbitrary effective range. Two different wave- 
functions are needed, depending on whether one is interested 
in the polaron [24, 35] or molecule [23, 36, 37] properties. 
The quasiparticle parameters for the polaron (energy E + and 
residue Z, effective mass) and the molecule properties 
can be found either variationally, or diagrammatically using 
the ladder approximation. Both approaches yield identical re- 
sults, which closely match independent Monte-Carlo calcula- 
tions [27]. The properties of the repulsive polaron, which is 
intrinsically unstable due to the presence of the molecule-hole 
continuum (MHC) and of the attractive polaron, are obtained 
from the self energy. In particular, the interaction induced en- 
ergy shift and the decay rate are given by the real part and 
twice the imaginary part of the self energy, respectively [12]. 

Previous treatments [12, 23, 24, 27, 35-37] were based on 
a universal scattering amplitude, describing broad Feshbach 
resonances. To include effects of the finite effective range we 
employ a many-body T-matrix given by [38, 39] 



T(K,co) = 



m, 



i -1 



2nh 2 a(K,co) 



-n(K,») 



(1) 



Here TiK = Pk + Pu is the total momentum with pu and Pk 
the momenta of 6 Li and 40 K , m r = mum^j (mu + '«k) the re- 
duced mass, Yl(K, co) the 6 Li - 40 K pair propagator in the pres- 
ence of the Fermi sea, and a(K, co) = «b g f 1 — b-b A e C m/<>Ii ) 
an energy-dependent length parameter, with a\, g , AB, Bq, 
and 8fl being the background scattering length, the width, 
the center, and the relative magnetic moment of the Fesh- 
bach resonance. £cm(K, co) = hco — h 2 K 2 / (2M) + Ef, with 
M = mu + ma, is the energy in the center of mass reference 
frame of the colliding pair. In vacuum and close to resonance, 
the scattering amplitude of our model has the usual low energy 
expansion 



fk 



- ik - r e k 2 /2 



(2) 



with the relative momentum /ik = (»?uPk — m KPu)/M. 
The effective range is approximated by r e ps — 2R*{\ — 
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ab s /a) 2 , where we introduce the range parameter R* = 
h 2 / (2m r a^hB8pL), see Ref. [22]. A detailed theoretical anal- 
ysis of this model will be given elsewhere [40]. 



2. Polaron peak in the spectral response 
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FIG. 5: Double Gauss analysis of the low-power spectra. The data 
are the same as presented in Fig. 2a plus additional data in an ex- 
tended range of \/{Kpd). (a) The Gauss function fitting the wide 
background is shaded grey. The Fourier-limited Gauss function, fit- 
ting the narrow peak, is coloured red (green) along the repulsive (at- 
tractive) polaron branch. We identify the narrow peak with one of 
the polaron branches only if its maximum signal exceeds a threshold 
value of 0.085, corresponding to two times the standard deviation of 
the noise in our data. Any smaller peak may be caused by fitting to 
a noise component. The lower panels show (b) the maximum signal 
of the narrow peak with the dashed line indicating the threshold, (c) 
the area under the wide Gauss function normalized to its maximum 
value, and (d) the detuning at the center of the narrow peak, provided 
that the peak signal exceeds the threshold, compared to the theoret- 
ical calculation of E + and £_ (red and green line). The error bars 
indicate the fit uncertainties. 



The spectra in Fig. 2a of the main text show a narrow, co- 
herent peak on top of a spectrally broad, incoherent back- 
ground. Here, we investigate these two spectral parts in more 
detail. Note that the background is actually better visible in 
Fig. 2b, but these spectra do not allow for a quantitative com- 
parison of the two parts because of the strong saturation of the 
narrow polaron peaks. 

The narrow peak stems from the attractive or repulsive po- 
larons, which correspond to well defined energy levels, pro- 
vided that the lifetime of the quasiparticle exceeds the pulse 
duration. As a consequence, the lineshape is expected to be 
Fourier limited except for the rapidly decaying repulsive po- 
larons very close to resonance. In contrast, the background 
is spectrally wide, on the order of Ep. The main contribution 
to the background stems from the MHC. Another contribution 
may arise from the excitation of additional particle-hole pairs 
in the Fermi sea when transferring to a quasiparticle with a 
momentum that is different from the momentum of the impu- 
rity in the initial state. 

We distinguish between the narrow peak and the wide back- 
ground by means of a double Gauss fit. Vertical cuts through 
Fig. 2a are presented in Fig. 5a together with the fit curves. 
The width a p of the Gauss function fitting the narrow peak is 
fixed to the one associated with the Gaussian fit of the Black- 
man pulse line shape used in the experiment, a p = 0.7kHz= 
0.019 Epjh. We constrain the width a), of the Gauss func- 
tion reproducing the background to 3 x 0.019 Bp/h < a\, < 
0.5 Ep/h. The lower bound avoids the misinterpretation of the 
narrow peak as background and the upper bound, correspond- 
ing to the maximal width of the continuum as obtained from 
the spectra in Fig. 2b, avoids unphysically large values of Ob 
when the background signal is weak. We find that the nar- 
row peak dominates for weak positive and negative interaction 
strength while the wide background dominates in the strongly 
interacting regime. This trend is shown in Fig. 5b and Fig. 5c, 
where we present the maximum signal of the narrow peak and 
the area of the background, respectively. Note, that the sig- 
nal in Fig. 5b is proportional to the area of the narrow peak 
since a p is kept constant. Figure 5d shows the detuning at the 
center of the narrow peak, which corresponds to the energy of 
the quasiparticles. The measured energies agree remarkably 
well with the calculation. The slight mismatch between the- 
ory and experiment may be attributed to systematic errors in 
the determination of Ep and Bq. 

The area of the wide background exhibits a maximum close 
to —1/ (Kpa) = 0, but it shows an asymmetry as it falls off sig- 
nificantly slower on the attractive (a < 0) side, see Fig. 5c. We 
attribute this asymmetry to the narrow character of the Fesh- 
bach resonance. The interaction becomes resonant when the 
real part of the inverse scattering amplitude, given in Eq. 2, 
is zero. This leads to the resonance condition = r e k /2, 
where a les is the value of the scattering length at which the in- 
teraction becomes resonant. In the limit of a broad resonance 
with r e = 0, this condition is fulfilled for any k at the center of 
the resonance, where the scattering length diverges. However, 
at a narrow resonance with r e < the condition requires a neg- 
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ative a res for k > 0. The mean square momentum in the Fermi 
sea is 3/5 x Kp, leading to a mean square relative momentum 
of 3/5 x (40/46 x Kp) 2 . Using this value for k 2 , and insert- 
ing r e as —2R* in the above resonance condition, we obtain 
— l/(fCffl) = 0.43. This represents an effective shift of the 
Feshbach resonance center, as we average over all momenta 
of the Fermi sea [41]. The magnitude of this shift agrees well 
with the observed asymmetry. Moreover, we find that many 
features at our narrow resonance appear to be shifted, e.g. the 
polaron-to-molecule crossing. However, the narrowness has 
many more implications and cannot simply be reduced to this 
shift. We will come back to this point in the context of the 
lifetime of the repulsive polaron, see Sec. 4. 



The repulsive polaron peak is clearly visible up to 
— 1/ (Kp a) ~ —0.3 while the attractive polaron peak vanishes 
already at — l/(Kpa) ps 0.9, see Fig. 5b. The fading out of 
the quasiparticle peak towards the strongly interacting regime 
approximately coincides with the position where the quasi- 
particle branches merge into the MHC. This shows that the 
polaron state is hardly observable as soon as it becomes de- 
generate with molecule-hole excitations. The MHC is not 
strictly limited to the range from E,„ to E„, — £f, as discussed 
in more detail in Sec. 3. It extends below E,„ — £f because 
of finite temperature effects. It also extends slightly above 
E m because of additional excitations in the spectral function 
of the molecules [13]. As a consequence, for finite tem- 
perature, the attractive polaron can become degenerate with 
molecule-hole excitations for values of the interaction param- 
eter above the calculated polaron-to-molecule crossing. This 
explains that the observed sharp peak is observed to disap- 
pear already at — l/(Kpa) w 0.9, which lies somewhat above 
the zero-temperature polaron-to-molecule crossing predicted 
at 0.6. 



It is interesting to consider the data analysis presented in 
Fig. 5b and c in relation to the common method of extract- 
ing the quasiparticle residue Z from the spectral weight of the 
narrow peak [42]. Close to resonance, we are in the linear re- 
sponse regime and our data can be interpreted in terms of this 
method. Our data suggests that this method leads to a signif- 
icant underestimation of Z. For example at —1/ (KfCi) ss 0.9, 
where the narrow peak of the attractive polaron vanishes, our 
theory still predicts Z rj 0.7. This underestimation is consis- 
tent with the one reported in Ref. [16], see also related discus- 
sion in Ref. [23]. A plausible explanation may be that such a 
method does not probe the polaron states alone, but also the 
molecule-hole excitations, which are degenerate with the po- 
laron state. Our alternative method of measuring the residue 
via the Rabi frequency, as presented in the main paper, offers 
the advantage of being much less affected by the molecule- 
hole contribution. In fact, only the coherent part of the quasi- 
particle is expected to produce Rabi oscillations, see Sec. 6. 



3. Molecule-hole continuum 

The spectra presented in Fig. 2b of the main text reveal the 
MHC. This continuum arises from processes where the rf field 
associates a 40 K impurity and a 6 Li atom out of the Fermi 
sea to a molecule. Here we present a simple model for the 
spectral line shape, which allows us to interpret the data up to 
— l/(Kp a) ?s — 1, see Fig. 6. 
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FIG. 6: Molecule association spectra for different values of the in- 
teraction parameter. The signal is the fraction of transferred atoms 
as a function of the rf detuning. The data correspond to vertical cuts 
through Fig. 2b. The dashed line is the line shape model for zero tem- 
perature and the solid line for finite temperature. The upper threshold 
of the theoretical spectra corresponds to E m . 

For modeling the line shape, we consider two-body pro- 
cesses in which the rf field associates one 40 K and one 6 Li 
atom to a molecule. Higher-order processes, involving more 
than two particles, are neglected in this model but are briefly 
discussed at the end of this section. Let us first consider the 
association of 6 Li and 40 K with momenta pu = p^ = 0. This 
results in a molecule at rest plus a Fermi sea with a hole in the 
center. The energy of this state is determined by the binding 
energy of the molecule and by the interaction of the molecule 
with the Fermi sea. It is given by E m and sets the onset of the 
MHC from the right (the top) in Fig. 6 (Fig. 2b). In general, 
6 Li and 40 K have finite initial relative momentum %k, leading 
to an initial relative kinetic energy in the center of mass frame 
E r = Ti 2 k 2 /2m r . The energy conservation of the association 
process is expressed in the Dirac 8 function in Eq. 3. As a con- 
sequence, the molecule spectrum extends downwards to ener- 
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gies below E m . We now consider an ensemble of 40 K and 6 Li 
atoms. Our experimental conditions are well approximated 
by a thermal cloud of 40 K in a homogeneous Fermi sea of 6 Li 
(see Methods). The momentum distribution of 6 Li is given by 
the Fermi-Dirac distribution f™(Eu), with Eu = Pu/2 m u- 
The one of 40 K is approximated by the Maxwell-Boltzmann 
distribution f^ B (E^), with E^ = p^/2m^. The latter distri- 
bution does not change its momentum dependence with po- 
sition, thus, no integration over space is needed to obtain the 
spectral response 



J?(A) « j j d 3 PLl d 3 p K jg>(Eu) 
f¥ B (E K ) J?(k) <?(-£,„+£, + A), 



(3) 



where &{k) is the Franck Condon overlap of the initial 
wavefunction with the molecule wavefunction. In our case 
the interaction in the initial state is negligible and &{k), 
as given in Ref. [43], reduces to &(k) « (E r /El) x l 2 (\ + 
E r /Et,)~ 2 . The parameter Ef, is the binding energy of a 
molecule in vacuum at a resonance with finite effective range 
and reads Ej, = h 2 / (2tn r a* 2 ) with the parameter [22] a* = 
— r e j (yl — 2r e /a — 1). In the calculation of JP(k), we do not 
account for interactions with the Fermi sea. Because of this 
approximation, we apply the model only for — 1 / (Kfd) < — 1 . 
For fitting the model line shapes to the experimental data, ad- 
justable parameters are the individual heights of the spectra 
and the center of the Feshbach resonance. The latter param- 
eter is required to be the same for all data sets in Fig. 6. In- 
dependently determined parameters are kgT /Bp = 0.16 and 
£f = hx 37 kHz. The model (solid lines) reproduces our data 
remarkably well. It allows us to pinpoint the resonance posi- 
tion to Bq = 154.7 19(2) G. This determination of Bo relies on 
our theoretical model to calculate E m . To test this model de- 
pendence, we replace E m simply by the binding energy of the 
molecule in vacuum plus the mean field energy, considering 
the corresponding atom-dimer scattering length [44]. Using 
this simple model, the fit yields a resonance position that is 
1 mG higher, which shows that the model dependence causes 
only a small systematic uncertainty. Moreover, the statistical 
fit uncertainty and the field calibration uncertainty are about 
1 mG each. 

For T — and all other parameters unchanged, the model 
provides the dashed lines in Fig. 6. The spectra show a sharp 
drop at A = E m — (40/46) E F , which corresponds to the asso- 
ciation of an impurity at rest and a majority atom at the Fermi 
edge. In an equal-mass mixture this process would occur at 
A = E m - (1/2) e F . Thus, the width of the MHC in the two- 
body approximation is much larger for a heavy impurity than 
it is for an equal-mass impurity and it is even narrower for a 
light impurity. 

The true zero temperature ground state is actually at the en- 
ergy E m — Ef, a molecule at rest formed from a 40 K atom at 
rest and a 6 Li atom at the Fermi edge. However, to reach this 
state, momentum conservation requires a higher-order pro- 
cess, i.e. the scattering of at least one additional 6 Li atom from 
and to the Fermi surface. Such processes are not included in 



the model presented here, which only considers the direct as- 
sociation of two atoms by an rf photon. 

In the strongly interacting regime the spectral function of 
the molecule shows additional excitations above the molecular 
ground state [13]. This leads to an extension of the MHC 
spectral response above E m , of which we find clear indications 
in our data. The lower panel in Fig. 6 shows finite signal above 
E m and the extension above E m is very evident in the strongly 
interacting regime, see Fig. 2b. 



4. Decay rate of the repulsive polarons 

We analyse the decay of the repulsive polarons by assuming 
that they decay into well defined attractive polarons or well- 
defined molecules. In this quasiparticle picture, the decay is 
associated with the formation of a particle-hole pair in the 
Fermi sea to take up the released energy. In this sense, the de- 
cay into the attractive polaron is a 2-body process and the de- 
cay into the molecule is a 3-body process. We calculate the de- 
cay rate for these two channels by including them into the po- 
laron self energy using a pole expansion of the 40 K propagator 
writing G(k,ffl) ~ Z+/ (Tico - E+ -h 2 k 2 /(2m K )) + Z_/(/zco - 
E- — h 2 k 2 / (2mic)) and a pole expansion of the T-matrix writ- 
ing T(k, t») ~ Z m g 2 /(h(o - (E m - e P ) - h 2 k 2 /{2M)). Here, 
Z± is the quasiparticle residue of the repulsive and attractive 
polaron respectively and Z m the quasiparticle residue of the 
molecule. The factor g 2 = 2nh A / (m 2 a* \ — 2r e /a) is the 
residue of the vacuum T-matrix for a general resonance. The 
details of this approach are given in Refs. [12, 45], the only 
difference being that here we include the effects of the finite 
effective range. The imaginary part of the self energy gives 
the decay rate of the wavefunction and we thus take twice the 
imaginary part to calculate the population decay. The 2-body 
decay into the attractive polaron and an additional particle- 
hole pair is calculated numerically to all orders in the T-matrix 
by inserting the pole expansion for the 40 K propagator in the 
self energy in the ladder approximation. For the 3-body de- 
cay into a molecule and an additional particle-hole pair, we 
include terms containing two 6 Li holes in the 40 K self energy 
[45], and an expansion to second order in the T-matrix relevant 
for —\/{Kpd) <C —1 yields 



pm ' 



tlKp 



64-Kpa Z\ 
45 7T 3 m^y/mu 

\ 5 



M 



(4) 



y/2(E+-E m + e F ) J a*yj\-2r e /c 



For simplicity, we have taken Z m = 1, which is an appropri- 
ate assumption for — 1 / (Kpa) -C — 1. The effect of the narrow 
resonance on the decay rate enters through the quasiparticle 
residue Z + , the energies E+, E m and directly through the 
effective range r e . This decay rate has the same a 6 dependence 
as the three-body decay in vacuum in the limit of a broad res- 
onance derived in Ref. [46]. The numerical prefactor however 
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FIG. 7: Decay rates of repulsive 4l) K polarons in a Fermi sea of 6 Li 
atoms, shown as a function of interaction strength (left) and of the 
energy of the repulsive polaron (right). Blue and red lines represent 
the two- and three-body contributions, respectively, while data points 
are the experimental findings as also shown in Fig. 3 of the main text. 
The results for the moderately narrow resonance under study here 
(solid lines) is compared with the theoretical results obtained for the 
universal limit of a very broad resonance (dashed lines). The experi- 
mental values of E+ are obtained by interpolation of the narrow peak 
position data Ap^, see Fig. 5d. 



differs since we have included the effects of the Fermi sea in 
a perturbative calculation. 

The results for the decay rates of repulsive polarons are 
shown in Fig. 7. The experimental data agree well with the 
theoretical results obtained for our narrow resonance (contin- 
uous lines) as already shown in Fig. 3 in the main text. For 
comparison, we also show the decay rates one would obtain in 
the limit of a broad resonance (dashed lines). We find that as 
magnitude of the effective range increases with respect to the 
interparticle spacing, the dominant two-body decay is strongly 
suppressed. This suppression is mainly due to a large reduc- 
tion of the attractive polaron residue Z . Instead, the weaker 
three-body decay increases, which we attribute to the reduc- 
tion of the polaron-molecule energy difference E + — E m + £f. 
Taking both decay rates together, the decay rate is at least an 
order of magnitude smaller at our narrow resonance as com- 
pared to the case of a broad resonance. It is important to note 
that this strong suppression of the decay at a given — \J{Kfcl) 
cannot be simply attributed to the effective resonance shift at 
our narrow Feshbach resonance as discussed in Sec. 2. When 
taking this shift into account, a suppression factor of five to 
ten remains. To highlight this point, we choose a representa- 
tion that is independent of the interaction parameter and that 
gives the dependence on the polaron energy, a direct manifes- 
tation of strong interactions. The right panel shows the same 
data and calculations as a function of E+. Also for a given 
E+, the repulsive polaron at our narrow resonance turns out 
to be much more stable than the repulsive polaron at a broad 
resonance. 



5. Decay of repulsive polarons to molecules 

The decay of the repulsive polarons, shown in Fig. 3 of 
the main text, is measured by applying a special three-pulse 
scheme (see Methods). In this section we exploit the flex- 
ibility of this scheme to study the decay to lower-lying en- 
ergy states in more detail. At a given interaction strength 

— 1/ (Kfd) = —0.9, we demonstrate that the repulsive polarons 
decay to molecules by showing that an rf spectrum taken after 
decay perfectly matches a reference spectrum of molecules. 

To populate the repulsive polaron branch, as done for the 
measurements of the decay rate, we tune the energy of the 
first pulse to E+, corresponding to A = 0.16£f at — \/{Kfo) = 
—0.9. The pulse duration (t p = 0.06 ms) and the intensity 
are set to correspond to a 7T-pulse in the noninteracting sys- 
tem. The second pulse removes the remaining non-transferred 
atoms by transferring them to a third spin state. In contrast to 
the decay measurement presented in the main text, we here use 
much more rf power for the third pulse to be able to efficiently 
dissociate molecules. For this purpose, we set t p = 0.3 ms and 
the pulse area corresponds to a 3 n -pulse in the noninteract- 
ing system. By varying the rf detuning, we record spectra for 
zero hold time (black squares) and for a hold time of 2 ms (red 
dots), see Fig. 8a. The peak at small positive detuning shows 
the back-transfer of repulsive polarons. The corresponding 
signal decreases with hold time, signalling the decay of the 
repulsive polaron. In addition, a wide continuum in a range 
of negative detunings rises with increasing hold time. Such a 
wide continuum involves coupling to high momentum states, 
signaling a short distance between 40 K and 6 Li. To confirm 
that this continuum stems from molecules, we compare it to a 
reference spectrum of the dissociation of molecules (blue di- 
amonds). We find a perfect match. To take such a reference 
spectrum, only the detuning of the first rf pulse is changed to 
directly associate molecules in the MHC instead of populating 
the repulsive polaron branch. We achieve a good association 
efficiency with A = — 0.54£f and t p = 0.5 ms. 

To study the evolution of the molecule population, which 
is fed by the decay of the repulsive polarons, we set the de- 
tuning of the third pulse to the peak of the molecule signal 
at A = — 13 Bp and record the signal as a function of the 
hold time, see Fig. 8b. A simple exponential fit yields a 
rate of about 1 ms~ 1 =0.0043 Ef/H, which is in good agree- 
ment with the measured decay rate of the repulsive polaron at 

— \/{kfo) = —0.9. The finite signal at zero hold time may 
have two origins. One contribution is some decay during the 
finite pulse durations of the three pulses, which are not in- 
cluded in the hold time. Another contribution may be the 
high momentum tail of the repulsive polarons as discussed in 
Ref. [16]. 

Note that we do not find any second sharp peak at negative 
detuning, which would indicate the population of the attrac- 
tive polaron branch. In case the repulsive polaron decays to 
the attractive polaron, the absence of the attractive polaron 
peak implies a very rapid subsequent decay of the attractive 
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FIG. 8: Decay of repulsive polarons to molecules at —\/{Kpd) = 
—0.9. (a) The black squares (red dots) show the spectrum right after 
(2 ms after) the repulsive polaron has been populated. The blue di- 
amonds show the dissociation spectrum of molecules for reference. 
The signal is the fraction of atoms transferred from the interacting 
spin state |1) to the noninteracting spin state |0). Note that the po- 
laron peak at positive detuning is highly saturated and thus its signal 
is not proportional to the number of polarons. (b) The rf energy de- 
tuning is fixed to A = — 1.3 Ef and the signal is recorded versus hold 
time. The error bars indicate the statistical uncertainties derived from 
at least three individual measurements. 



polaron to the MHC. Such a fast decay of the attractive po- 
laron to the MHC is consistent with the very small signal of 
the attractive polaron peak throughout the regime of strong 
interaction as discussed in Sec. 2. 

Let us briefly discuss the possible role of inelastic two-body 
relaxation in the 6 Li - 40 K mixture, which is energetically pos- 
sible as 40 K is not in the lowest spin state. This process was 
identified in Ref. [21] as a source of losses. However, this re- 
laxation is about an order of magnitude slower than the mea- 
sured decay rate of the repulsive polaron and thus does not 
affect our measurements. 



6. Rabi oscillations and polaron quasiparticle residue 




FIG. 9: Linear increase of the Rabi frequency CI with the unperturbed 
Rabi frequency The left (right) panel shows the driving to the 
repulsive (attractive) polaron. The solid lines are linear fits without 
offset and demonstrate the proportionality Q. °c £2q. 



pends on the matrix element of the rf probe between the ini- 
tial state |0) and the final state |1). Since the probe is ho- 
mogenous in space, it does not change the spatial part of the 
atomic wavefunction and it can be described by the opera- 
tor [47] R °c Q. Q £ q (fl| q fl 0q + h.c.) where d J q (a, q ) creates (an- 
nihilates) a 40 K atom with momentum q in spin state i and 
Q.q is the unperturbed Rabi frequency of the |0) to |1) tran- 
sition in the non-interacting case. Considering for simplicity 
an impurity at rest, the initial non-interacting state is given by 
\I) = flg q=0 |FS) where |FS) is the 6 Li Fermi sea. The final 
polaronic state at zero momentum can be written as [35] 



VZa\ q=0 \FS) 



E < 

q<hK F <p 



j<V P £pNIFS) 



(5) 



where b q (b q ) creates (annihilates) a Li atom with momen- 
tum q. The second term contains a Fermi sea with at least 
one particle-hole excitation and thus is orthogonal to an un- 
perturbed Fermi sea. Therefore the matrix element reduces to 
(F\R\I) = \/ZOo and we obtain the Rabi frequency 



(6) 



We neglect the momentum dependence of the quasiparticle 
residue and do not perform a thermal average over the initial 
states, which we expect to be a good approximation since T <C 

£f Ag- 
in Fig. 9 we plot the observed Rabi frequency £2 as a func- 
tion of the unperturbed Rabi frequency £2o. We find that 
the proportionality £2 « £2o holds over a wide range of rf 
power. The measurements presented in the main text, taken at 
£2() = 2n x 6.5kHz and 12.6kHz, are safely within this range. 



For high rf power, the signal is well beyond linear response 
and the 40 K atoms exhibit coherent Rabi oscillations between 
the spin states |0) and |1). In this regime the oscillations are 
so fast, that the polaron decay plays a minor role and can 
be ignored to a first approximation. The Rabi frequency de- 
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